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Abstract 

In the theory of deep inelastic scattering (DIS) the final state interaction (FSI) 
between the struck quark and the remnants of the target is usually assumed to be 
negligible in the Bjorken limit. This assumption, still awaiting a full validation within 
nonperturbative QCD, is investigated in a model composed by two relativistic particles, 
interacting via a relativistic harmonic oscillator potential, within light-cone hamilto- 
nian dynamics. An electromagnetic current operator whose matrix elements behave 
properly under Poincare transformations is adopted. It is shown that: i) the parton 
model is recovered, once the standard parton model assumptions are adopted; and ii) 
when relativistic, interacting eigenfunctions are exactly taken into account for both the 
initial and final states, the values of the structure functions, averaged over small, but 
finite intervals of the Bjorken variable x, coincide with the results of the parton model 
in the Bjorken limit. 
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1 Introduction 



The parton model proposed by Bjorken and Feynman almost thirty years ago |]I|, |^ has 
turned out to be a convenient language for discussing different properties of deep inelastic 
scattering (DIS), although, according to the theory based on the operator product expansion 
[P, the parton model, even in the Bjorken limit, is accurate only up to anomalous dimensions 
and perturbative QCD corrections. 

If q is the momentum transfer in DIS, Q = Iq'^I^^^, P' is the four- momentum of the 
target in the initial state and m is its mass, then the Bjorken limit is defined by the condition 
that Q ^ m and the Bjorken variable x = Q'^/{2P'q) is not too close to and 1. In the 
parton model, it is assumed that the final state interaction (FSI) of the struck quark with 
the remnants of the target is a higher twist effect, i.e., an effect which is suppressed at least 
as m? /Q"^. The qualitative motivation of this assumption is that the time needed for the 
absorption of the virtual photon by the struck quark is much smaller than the time of its 
hadronization and therefore in the process of absorption the struck quark can be considered 
as approximately free. It has been also shown within the framework of the coUinear expansion 
1^ that the Feynman diagrams describing the FSI are indeed suppressed at high Q^, once 
proper assumptions are adopted. 

It could be thought very surprising that FSI can be disregarded in presence of confine- 
ment (see, e.g., Ref. [Q). Indeed, at first glance the picture which follows from confinement 
fully differs from that given by the parton model, i.e., no interaction in the final state and 
a continuous spectrum. In particular, while the structure functions in the parton model are 
continuous, in models where confinement is taken into account the structure functions are 
linear combinations of delta functions. However the two models can be reconciled within the 
framework of the theory of distributions. An analogous situation takes place in the case of 
the reaction e^e~ hadrons (see, e.g., the discussion in Ref. Q and references therein). 

Of course the relevance of the FSI has to be studied in the framework of nonpertur- 
bative QCD, but in absence of a full solution it is desirable to consider models in which the 
structure functions can be calculated exactly and therefore it is possible to check whether 
the FSI is indeed a higher twist effect. 

In the last years several nonrelativistic two-body models of inelastic scattering were 
considered (see, e.g., Refs. [0, §]) and it was shown that the effects of the FSI in these 
models are in agreement with the standard parton model assumptions, even in the case of 
confining interactions. However, DIS in the Bjorken limit by no means can be considered 
nonrelativistically. Indeed, in this limit the energy transfer g° in the target rest frame has 
to satisfy the condition g° ^ m, while the nonrelativistic approach holds only if g° <^ m. 

The FSI for a relativistic particle bound in an external field was considered in Ref. 
0. However, in this case one has a two-body problem, where one of the particles (the source 
of the field) has an infinitely large mass. Therefore this problem cannot be studied in the 
Bjorken limit, where Q is much larger than all the masses involved in the problem. Another 
relativistic approach to the FSI has been considered in Ref. in the framework of the 
Bethe-Salpeter formalism, but it was assumed, as usual, that only Feynman diagrams from 
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a certain class (so called handbag diagrams) are dominant. Then the result of Ref. 0] that 
the FSI is a higher twist effect has been confirmed. 

The aim of this paper is to investigate the role played by FSI in DIS for an exactly 
solvable relativistic constituent quark model (CQM), within the light-cone (= front-form) 
hamiltonian dynamics. As well known, in the framework of relativistic CQMs (see, e.g., 
Ref. [§]) confinement can be ensured by choosing a quark-quark potential such that the 
mass operator of a system with a fixed number of relativistic constituent quarks has only 
the discrete spectrum (while in QCD confinement is understood as the property of the quark 
and gluon Green functions to have no poles for real values of the mass). Our purpose is to 
verify whether the naive treatment of confinement in relativistic CQMs is compatible with 
the parton model. 

We consider a simple system composed by two relativistic particles interacting via 
the relativistic harmonic oscillator potential. We adopt an electromagnetic current whose 
matrix elements exhibit the correct properties under Poincare transformations and fulfill the 
current conservation, as shown in In the proper Breit frame, the relevant components 



of the current are the same as in the parton model. Then, in the framework of the light- 
cone hamiltonian dynamics, we can derive exact expressions for the DIS structure functions, 
including the FSI effects, and show that in the Bjorken limit the exact results coincide with 
those given by the parton model, after an average over small intervals of the scaling variable 
X has been performed. This average features the finite detector resolution and allows us to 
avoid some mathematical technicalities of the theory of distributions. 

The paper is organized as follows. In Sec. ^ we explicitly define a model describing 
the interaction of two relativistic spin 1/2 particles, then in Sec. ^ we calculate the DIS 
structure functions of this system with the standard parton model assumptions and in Sec. ^ 
the structure functions are calculated using the exact two-body wave functions. Conclusions 
are drawn in Sec. 5. 

2 Relativistic harmonic oscillator potential 

We consider a system of two different particles with the same mass, mo, and spin 1/2. 
To describe such a system it is necessary to choose first an explicit form of the unitary 
irreducible representation (UIR) of the Poincare group pertaining to each particle. There 
are many equivalent ways to construct an explicit realization of such a representation [ p^ . 
We choose the realization in the front form of dynamics (see, e.g., Refs. [|lli |l^, |13|, |l^, |15 



Let p be the particle 4-momentum, s be the spin operator and a be the spin projection 
on the z axis {a = ±1/2). We define = (p° ±p')/V2, and we use p± to denote the 
projection of p onto the plane xy. The one particle Hilbert space can be chosen as the space 
of functions (p{p±,p'^, cr) such that 
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The Hilbert space Ti for the representation of the Poincare group describing a system 
of two free or interacting particles is reahzed in the space of functions Pi^, cii, p2±, vt i ^2) 

such that 

0-10-2 i=l 

Instead of the variables pi^, pf, p2±, pt, let us consider the variables P±, P~^, and k, 
where Pi = pi± + p2±, = Pi + pt s-^d the relative momentum k is defined as follows 
T^ , |T3[. We first define the quantities 

^=^, k±=Pi±-^P± (3) 



and then define k = {k±,kz), where 



^_1 = _^, uj{k) = {ml + kY\ (4) 



A direct calculation shows that 



Y[dpi{p., 

i=l 

dpiint) 



zj.,pt) = dp{P±,P^)dp{int), 
dk±d^ dk 



2(27r)3e(l-0 (27r)3cu(^) 



(5) 



If the particles do not interact with each other, the generators of the two-particle 
representation are equal to sums of the corresponding one-particle generators fl^. The 
result is that the free-mass operator of the system is Mq = 2uj{k) and the two-body spin 
operator is equal to 

S = U~\l{k) + si + S2]U, U = v{k, si)v{-k, S2), (6) 

where v{k,s) is the Melosh matrix [|T6|. In the given context the Melosh matrix was first 
considered by Terent'ev fl^ : 



^ ^ 2ie is^k^ 
v{k,s) = exp{ — -J arctg 



mo + uj{k) + kz 



(7) 



with fc_L = |fc_L|, j = 1, 2, / = 1, 2, and ei2 = -621 = 1, en = €22 = 0. 

Let us define the "internal" Hilbert space Tiint as the space of functions x = x{k, cti, (72) 
such that 
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J2 J \xik,ai,a2)\'^dp{int) < oo 



If particles 1 and 2 interact with each other, then the representation space is the same, 
but the generators differ from the free ones (see, e.g., Eqs. (A34-A36) of Ref. [|13| or Eq. 



(2.34) of Ref. [|T5l). The interaction should be introduced in such a way that the new set 
of generators, as well as the set of free generators, satisfies the commutation relations of 
the Poincare group Lie algebra. By analogy with Ref. [1^, this can be done by replacing 



the free mass operator Mq by an interaction dependent mass operator M, which acts only 
through k, ai, 02 and commutes with the spin operator S given by Eq. (^. As well known, 
the system of generators obtained in such a way defines the front form of dynamics |]TB|, in 



which only three generators are interaction dependent and the other seven generators are 
free. 

Then, let us define the internal space lUinf, whose elements \l/(/c, ai, a^) = {k, ai, cr2|\E') 
are normalized as in the nonrelativistic quantum mechanics: 

5: /|v^(fc,ai,a2)p-^ <oo (9) 
0-10-2 •' v^^/ 

As follows from Eqs. (|^), (|]), (^ and @, if the relation between the spaces Hint and 7i-„j 
is defined as 

Xik,a^,a2) = {k,a^,a2\x) = {k, a^, a2\U-'\^MkY/^ (10) 

then the two-body spin operator in Ti'^^t has the standard "non-relativistic" form Snr = 
ilk) + si + S2 (that also coincides with the instant-form one, see, e.g., |T^, |T3|). 

Now we use the observation which is the essence of the " minimal relativity principle" 
T^■. if M is the mass operator in 7i-„^ and the interaction operator V in TC'^^^ is defined as 



= Mq + V, then the equation M^\I^„ = M'^'^^ for the eigenvalues M„ and eigenfunctions 
of the operator M has the same form as the nonrelativistic Schroedinger equation in 
momentum representation: 

P ^ 

( — + V)^n{k, ai, (T2) = K^„(A;, ai, as) (11) 
mo 

where 



V = V/Amo, En = {Ml - 4m2)/4mo (12) 

The operator V should satisfy the same conditions as in the nonrelativistic quantum me- 
chanics: the operator (fc^/mo) + V should be selfadjoint and V should commute with Snr- 

We can formally introduce the operator f = id/dk which is canonically conjugated 
with k. It is well known that in the relativistic case there is no operator which has all the 
properties of the position operator. Therefore r has not all the properties of the operator 
of the relative radius-vector between particles 1 and 2; it has such properties only in the 
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nonrelativistic or classical limit. Nevertheless we can choose the operator V as the operator 
of multiplication by a function V(r), where r = |r|. In particular the function V(r) can be 
chosen in such a way that the operator M (and hence M) has only a discrete spectrum, i.e., 
the property which in CQMs is associated with confinement. 

We choose the function V(r) in the form V(r) = a^r^/mo, where a is some constant 
with the dimension GeV . Then Eq. (|Tl|) is the well-known equation for the harmonic 
oscillator. We stress, once more, that in Eq. (0) the mass operator is a fully relativistic 
one, though it coincides with the Schroedinger operator. The solutions of Eq. ([TI|) are given 

by the products ^n,a[,u'2 = i^n{k)ipa[,a'^{<yi,(^2), where ipa[,a'^{<yi,cr2) is the spin eigenfunction 
of particles 1 and 2 with eigenvalues (T[,a2, respectively, and 

a a 
1 

i = x,y,z. (13) 

In Eq. (|13D rij = 0,1,2..., and H^n is the Hermite polynomial of the rij-th order. The 
normalization of ipm is S'^oo\'4^n,{ki)? = 1- As follows from Eqs. (|ll]) and (|l|), the 
eigenvalues of the mass operator are equal to 

Mn = 2[ml + a^{2n + ?,)Y'^ (14) 

with n = + Uy + Uz- 

In order to calculate the DIS structure functions we have to know the wave function 
of the two-body system in the space Ti when the internal wave function Xn,u[,(j'2 in the 
space l-iint is related to '^H,cT[,a'^ according to Eq. (0) and the system as a whole is in the 
eigenstate of the operators P± and P+ with the eigenvalues P' ^ and P'+, respectively. Let 
us first normalize as follows the one-particle eigenstates with the four-momentum p' and spin 
projection a'\ 

= 2{2Txfp+5^^\pl - p\)5{p + - p+)5„,^^, (15) 

(^ctV is the usual Kronecker symbol). Then the state |-P'±, -P'"'', Xhcj' a ) described by the 
wave function 

P'+5(2)(p"^ _ P'^)<5(P"+ - P'+)x„v;,.^(^,ai,a2) (16) 

From Eqs. (p!0|) and (p!3| ) it is clear that | |Xn,o-i,o-2 1 1 = 1; if the spin eigenfunctions are 
properly normalized. 
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3 DIS structure functions in the parton model 



Let us consider in some detail how the structure functions can be obtained within the 
parton model in our light-cone framework. The results are well known (see, e.g., Ref. ||20||), 
but the purpose of the following derivation is to make a comparison with the one given in 
Sec. 4, where the final states are described by the exact harmonic oscillator wave functions. 
For simplicity we will consider the case where the total spin of the initial system is S" = 
(this is by no way a restriction: indeed in the case where S = 1 one can obtain the same 
results as well). Then the internal wave function of the initial state will be 

— * — * 

Xo{k,ai,a2) = {k,ai,a2\xo) 

= uikY/^Mk){ai,a2\U-'\ips=o) (17) 

Let J^{x) be the electromagnetic current operator for the system under consideration 
(/i = 0, 1,2,3), where x is a point in Minkowski space. It is well known that if the initial 
state is |-P'±, -P'"*", Xo) then the DIS hadronic tensor is given by 



W^" = — f exphqx) ■ 
Atx J 

■{P'±, P'^, Xo\J''{x)rmP'±, P'+, Xo)d'x 



The coordinate dependence of the current operator is fully defined by translational invari- 
ance, according to which 

j/^(a;) = exp{iPx)J^'{0)exp{-tPx) (19) 

where P is the total four-momentum operator. Therefore, as follows from Eqs. ( [T8|) and 
(0), 

47r ^ 

(P'x,P'+,Xo|J''(0)|X)(X|r(0)|P'x,P'+,Xo) (20) 

where a sum is taken over all possible final states \X) and Px is the four-momentum of the 
state \X). It is also well known that, as a consequence of Poincare invariance and current 
conservation, the unpolarized hadronic tensor has the form 

W^''{P',q) = - gnPi{x,Q) + 

_L.(p._^,(P'-^)^.(..«). (21) 
where g^'^ is the Minkowski tensor. 
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In the parton model one assumes that in Bjorken hmit the operator J'^(O) can be 
taken in the impulse approximation, i.e., J^(0) = (0) where the current operator 

Jf (0) acts as follows 

where are the Dirac 7 matrices and w{pi, (Ti) is the light-cone Dirac spinor. 
The following representation for the 7 matrices has been adopted 



(22) 



7 






1 




1 




-Ti 


1 





, 7^ = 


-1 


Ti 



(23) 



where i = 1,2,3 and Xj are the Pauli matrices. Then the light-cone Dirac spinor can be 
written as 



(24) 



where xi'^) is the ordinary spinor describing the state with the spin projection on the z axis 
equal to a and the matrix P{g±,g~^) (where g = p/rrio) has the components 

Pu = P22=2'^\9^Y^\ /5i2 = 0, [321 = {g. + igy) P22 (25) 

For simplicity we assume that particle 1 has unit electric charge and particle 2 is 
chargeless. With this assumption, from Eqs. (||), (|l^) and (|22D , we have 



(26) 



where ^' and k' are constructed by means of Eqs. (0,§) from the vectors p[ and P', where 



P2 



P'lL = P'±-P2±,Pl =P^ 

As noted above, one of the major parton-model assumptions is that the FSI of the 
struck quark with the target remnants can be neglected. This implies that in our case the 
final states X in Eq. (BO) are the states of two free particles: 



|X) = bi",^i")b2",a2") 
and therefore Eq. ( pOD for the hadronic tensor can be written as 

^''^ = ^ E /(2vr)V^)(P' + g-P")- 

o"l"a-2" 

(P'±,P'+,Xo|^'^(0)b;ai",p;,a2") 

(p;,ai",p2,a2"|r(0)|P'^,P'+,Xo) 
2 

Y[dpi{pl^,p+) 

i=l 



(27) 



(2J 
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where P" = p'^ + V^- turn, as follows from Eqs. (|), (0), (^) and (|27|) , Eq. (pSD can be 
written as 

where ^ = ^^^7^ and, as follows from Eqs. and (^), 

= (i-e')p'+ = (i-r)(p'+ + 9+) (30) 

Therefore the relation between {k'j_,^') and (^1,^") is 

fcl = + ie - c)p'± + (1 - cip'^ + 9+) = + (si) 

In the parton model the initial hadron is considered in the infinite momentum frame 
(IMF), i.e., in a reference frame where P'^ is positive and very large. It is also assumed that 
the transverse momenta of the constituents in the initial and final states are restricted by 
some value (say 300MeV/c). In view of this, a suitable choice of the reference frame is such 
that P'± = q± = 0. Then, from Eqs. (^) and (|3TD , one has 



Pi± = P'i± = kl = k'^, = -k'± (32) 

The above conditions do not define the reference frame uniquely, since one can still 
boost this frame along the z axis and choose the Breit frame, where P" + P' = 0, in order to 
study DIS (let us recall that such boosts along the z axis are kinematical in the front form). 
It is important to point out that, as shown in Ref. in the Breit frame the one-body 

current operator, Eq. (pB]), is fully compatible with the Poincare trasformation properties. 

In the chosen frame, as follows from the definition of the Bjorken variable x, one has 
in the Bjorken limit 



gO = 2P;(1 - x), P'+ = V2P;, q+ = -V2P'^x, 

P"+ = V2P^{l-x) (33) 

and the relation between ^' and ^" (see Eq. (|3TD) becomes 

r = x + (l-a;)r (34) 



9 



Furthermore, by a direct calculation using Eq. (p3D, one obtains that the relevant ele- 
ments of the hadronic tensor (PT|) are given in the Bjorken limit by the following expressions 



W^' = 5^'F,{x,Q), {j,l = 1,2) 

= -^J—^[F2{x,Q) - 2xFi(x,Q)], 

W+- = W-+ = ^[F2(x,Q) -2xFi(x,Q)], 

W-- = '^[F2{x,Q)-2xF^{x,Q)] (35) 



Finally, in the Bjorken limit the argument of the delta-function in Eq. (pQ]) becomes 
proportional to 

Q'il-x) (ml + kl^) 



X r(i-r) 

Therefore we have two solutions for the first of which is close to and the second is 
close to 1. A direct calculation based on Eqs. p^), (p5|), (p9D, (|32D and (Q) shows that in 
the Bjorken limit the contribution of the second solution is negligible. Indeed, when ^" is 
close to 1, then ^' is close to one too (see Eq. (0)) and k'^ goes to infinity (see Eq. (|)), 
making vanishing Xo iii Eq. (^91) . In what follows we describe only the contribution of the 
first solution (^" close to 0). 

The only large component of the momentum of particle 1 in the initial state is p'l' and 
the only large component of the momentum of this particle in the final state is pi~ . Then, 
as easily seen from Eq. (^), the spin structure of the tensor W^" in the Bjorken limit is 

~ Tr{-f-Yl^Y) (36) 

where Tr stands for trace. This implies that all the longitudinal components of W'^'^ are 
equal to zero, i.e., 

W^" = W-" = = W"- = 0, (37) 
since (7"^)^ = (7")^ = 0. The only non-vanishing components are the transverse ones, viz. 

W^' = /xo(^l,a:,ai",a2")*M.,.vi+2ze,,(s^).iX]- 

a-i"cr[a2" 

Hp 

In obtaining Eq. (|38| ) we have used the fact that, if ^" is close to 0, then in the Bjorken limit 
^' = X (see Eq. (0)). As follows from Eq. (|^), the quantity ^' is the momentum fraction 
of particle 1 in the initial state. Therefore we obtain the well-known result that the Bjorken 
variable x has the meaning of the momentum fraction of the struck quark in the IMF. 
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The structure function Fi{x,Q) is given by the symmetrical part of W^'' and the 
structure function F2{x,Q) by the longitudinal components. Let us introduce the notation 



dk I 



0"l,0-2 



2(27r)3x(l - x) 



dk 



2(27r)3x(l -x) 



(39) 



where (see Eqs. (0) and (|ll) Xo(A;±,0 = uj{kf'^i)o{k). Then, as follows from Eqs. (D, (|) 
and (H), p{x)dx is the probability to have a momentum fraction of particle 1, in the initial 
state, falling in the interval {x,x + dx). By using Eqs. (|35D , (W^), (|38| ) and (|39|) , one obtains 



(40) 



and the Callan-Gross relation [22 



F2(x) = 2xFi(x) 



(41) 



Therefore in the Bjorken limit the structure functions Fi and F2 do not depend on Q, namely 
one has Bjorken Scaling. 

4 DIS structure functions with the relativistic har- 
monic oscillator wave functions 

Let us consider the exact hadronic tensor for two particles, interacting via the rela- 
tivistic harmonic oscillator potential introduced in Sec. 2, and let the initial state be the 
5* = state as in Sec. 3. Therefore in Eq. (|20|) we use, as final states |X), the exact 
eigenf unctions defined by Eqs. {^), (|13]), (|16]). Then the exact expression for the hadronic 
tensor of the system under consideration is 



{P'^, P'+, xol J'^(0)|P;, P" + , XnVi",-2") • 

dp(p;,p"+) (42) 

where the four- vectors P"„ = P' + q have the components P"^ = P""*", P" n± = P±, 

P"- = (M2 + p;2)/2P"+, 

|XnVl",<.2"(^)) = Cu(^)l/Vn(^)f/~V.i",.2"), (43) 
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and I ,0-2) is the normalized spin eigenstate of particles 1 and 2 with the eigenvalues cJi, (T2, 
respectively: 

llv^^xi.aJI = 1. (44) 

We will calculate the tensor ( ^21) in the same reference frame as in the preceding 
section, i.e., in the Breit frame (P" +P' = 0) with P'± = q^ = 0. We assume J^{0) = Jf (0) 
for /i = +,1,2 and use current conservation to define Ji{0). As it is shown in |^T|, this 
current operator is an allowed choice, which is fully compatible in our reference frame with 
Poincare transformation properties, and trivially fulfills current conservation. Then, from 
Eqs. (P, (i, ®, (|T]) and (||), 



cri"cr2" nxTiyUz 



Ml 



(Ti crj 



2(27r)3^'f'(l-^" 



(45) 



where m = 2[mQ+3a^]^/^ is the ground state eigenvalue of the mass operator and the relations 
between the quantities p'i,p"i, k'j_,^', k'^,C' are given by Eqs. ( pOj ) and (PTf). In the reference 
frame under consideration, k'l_ = k'j_ as in the parton model, and the relation between ^" 
and ^' is still given by Eq. (p^). 

The spin sums in Eq. (|45|) can be easily done taking into account the following 
expressions for the Clebsh-Gordan coefficient in the initial state, the Melosh rotation and 
the matrix element of the parton current, viz. 



(aiaalOO) = i{ai\-^\a2) 

—* — * 

_^ _ mp + Lo{k) + k^- i2s- (e^ x k±) 
^2{u{k)+mo){u;{k) + k,) 

(o-il q^iicTym + kr," + iky'cr^) + 2p'{^k^" \a[) 

Pi"Vi+ 



(46) 



Then Eq. (^) becomes 



Xn(^±,0 



1 



mo + 



k^ 



uj{k) + mo 



dk±d^ 



2(27r)3e^(l-0 
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1 








1 







[u;(A;) + mo] 
q^kykx 



2(27r)3e'e(l-0 
dk±d^ 



2(27r)3ee(l-0 

dkj_d^ 1 2 
2(27r)3e'e(l-0 



+ 



where in order to simplify the notation, the variables /c^ and ^" of Eq. 
replaced by k± and ^; moreover (see Eqs. (0) and (^31) ) 



(47) 



have been 



Xnik) = uoikf'^i^S), Xoik') = uj{k'f'^^o{k'). i' = x + {I - x)i, 
and the relation between k {k') and {k±,^) {{k'j_,C,')) is given by Eq. (^). 
4.1 A pedagogical example 



(4J 



In order to begin with a simple, pedagogical model that allows us to introduce the math- 
ematical tools for the general case (to be considered below), we assume for the moment 
a <^ tuq. This assumption implies that the relevant momenta k' in the wave function of 
the ground state satisfy the condition \k'\ ^ rriQ (see Eq. (pTsP), but does not destroy the 
relativistic nature of the highly excited final states. Then Eq. (^) greatly simplifies, putting 
in evidence the relevant integral over ^. With the help of Eq. (pSf), Fi{x, Q) becomes 



F^ix,Q) 



2^2 



mQX 



2(1 -a;), 
Xn(A-±,Oxo(A;±,^' 



J2 Sirn^ + 



X] 



X 



dk±dC, 



(49) 



So far we have not used the explicit expressions for the harmonic oscillator wave 
functions and therefore, if the FSI is neglected in Eq. (^9]), the result for Fi{x, Q) should be 
the same as in Eq. (|40|) . Disregarding the FSI in Eq. (^9|) implies the following replacements 



^ r dkldC ^2 ml + k"f 

7 2(27r)3r(l-r)' " ^(1-^")' 

Xn{k±,0 ^ 2{2nfC{l-C)S^'\k^-kl)6{^-n (50) 



Then it is easy to see that we again arrive at Eq. (|40|) for Fi{x, Q). 

Let us go back to Eq. (|49|) , where the exact final state wave functions are used. Since 
only the values ^ mo are important in Eq. (^Qj), we can neglect k± in u^k). Then, as 
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follows from Eqs. @, (pTsD and (^81) , the integration over k± in Eq. (^91) is trivial and as a 
result 



2^2 



\ml + feD^/Vn(A:.)xo(n ,^(^.^)3f2(i_^) r (51) 

where n replaces ra^, Xo('C') = ^^e relation between ^ and is given by Eq. (^) 

and one-to-one relations C, ^ kz and ^' /c!. can be obtained from Eq. @. In particular if 
the dependence of uj{k) on k± is neglected, then the relation between ^ and A;^ is given by 

^1 K , mo(e-l/2) 



g2(l-x) 



2^(1-0 K + A;2)i/2 

and the relation between ^' and A;^ is the same. 

As a consequence of Eq. (|5^) , Eq. (|5lD can be written in the form 



(52) 



2^2 



|2 



(53) 



This is the last stage where we still can return back to the parton model. Indeed, as 
follows from Eqs. (pOD and (|5lD, neglecting FSI in Eq. (^31) implies the following replacements 



<;v(t) ^ 2T(mJ + kiy"d(h - K) (54) 



With these replacements it is easy to see that, in the Bjorken limit, Eq. (|53D leads to Eq. 
(|40| ) for Fi{x, Q), if Xo(^) is given by Eq. (148|) and a ^ mo. 

In the last part of this section we will consider Eq. (|5^) with the function ipnikz) 
given by the exact eigenfunction, i.e. by Eq. (pTSl). 

As follows from Eqs. ( [T4|) and (^), n = Q^(l — x)/8a^x and therefore, if Q is large, 
only large values of n are important in Eq. (^3[). Taking into account Eq. (0) we can write 
Eq. (|53D in the form 
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where 



f(n x) - f-ir r (5Q) 

f[n,X)-[ L) y_^^,(^,^)i/2^^2 + ^2)1/4 ^^^) 

We will show in Appendix A that f{n, x) is a smooth function of x and has a finite limit for 
n oc. 

Now the following question arises. While in the parton model Fi{x) and ^2(0;) are 
continuous functions of x, it is clear from Eq. (^) that at fixed Q the function Fi{x,Q) is 
a linear combination of delta-functions, which are not equal to zero only for discrete values 
of X. As noted in Sec. |l], the first impression is that a correspondence between discrete 
and continuous cases cannot exist. However, Eq. (|55|) is meaningful only in the realm of 
the distributions. Within such a framework, the correspondence between the discrete and 
continuous cases could be shown in the Bjorken limit, since the discrete values of x, where 
Fi{x,Q) is not zero, become closer and closer as Q increases. 

In order to simplify the mathematical discussion, let us note that experiments allow 
one to determine not the very function Fi{x, Q), but its average values over some bins in x 
and Q. Therefore, let us consider the integral 

Q) = — i — r Q)dx, (57) 

X2 — Xi Jxi 

where x belongs to the small interval [xi,X2], such that X2 — xi ^ x. It is clear that for 
large Q there exist many values of n such that Q'^/{Q'^ + 8a^n) G [xi,X2], even for a small, 
but finite value of X2 — Xi. Therefore the integral (pT]) is a smooth function of Q. Let us 
define ni and n2 as follows 

ni = — , n2 = — — . 58 

Sa'^Xi Qa'^X2 

As a consequence rii — n2 = {x2 — X\)n2l\x\{\ — 0:2)] ^ n2. From Eqs. (^) and ( |^ one has 

87r^(l - x){x2 - Xi)Q^ 

The scaling property of the function Fi holds if f{n,x) —>■ g{x) when n 00. As a 
matter of fact, in this case the quantity SnL„2 /('^^ large values of rii and n2 becomes 

{rii — n2)g{xY and therefore from Eq. (|58D one has 



Moreover, as follows from Eqs. (|39D , ( ^OD and (|5^) , the result (pOf ) for the structure function 
will be equal to the parton model one if 

hm /(n, x) = g{x) = ry^ (61) 
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where Xo{^) = fT\>^'^'^o{kz{x)), with kz{x) = mo(a; — l/2)[x{l — x)]~'^^'^. Indeed in our case, 
where the initial state is a harmonic oscillator ground state with a <^ mo, the structure 
function Fi{x) in Eq. (|40| ) becomes 

Fiix) = (62) 
87rx(l — X) 

In Appendix A it is shown how Eq. (plf ) can be proved. In conclusion our result for 
Fi{x,Q) (Eq. (^5D) is indeed compatible with the parton model, once an average over bins 
of X is performed. 

As follows from the analysis of Appendix A (see in particular the discussion about 
f2{n,x)), the main contribution to the hadronic tensor is given by the region where kz is 
negative and \kz\ is very large (recall that t = kz/a). In this region ^ is small and — in the 
Bjorken limit. In turn, the quantity ^' is close to x (see Eq. (|^)). We see that the usual 
interpretation of the Bjorken variable x as the momentum fraction of the struck quark is 
still valid if the FSI is not neglected. 

We recall that, with our choice of the reference frame, and \pij_ \ can be neglected 
with respect to mo (see Eq. (^O])) in the evaluation of the matrix elements of between 
light-cone Dirac spinors. Therefore a direct calculation using Eqs. ( p^ and ( |25D shows that 
the matrix elements w{p[,ai")'y~^w{p[,a[), up to constant quantities, are proportional to 
^(ti"<t[{p"~^P~^Y^'^ ■ Then, as shown in Appendix C, W~^~^ is vanishing in the Bjorken limit, 
and from Eq. (|35|) one obtains again that the Callan-Gross relation (^11) holds. 

Finally, using once more Eq. (|55D, one obtains that all the longitudinal components 
of the hadronic tensor become negligible in the Bjorken limit, with the exact wave functions 
for the final states as well as in the parton model. 



4.2 The general case 



In the general case, i.e. for any value of the harmonic oscillator strenght a, the parton model 
can be recovered once again. The starting point is Eq. (^) that is exact. After averaging 
over the x-bins, as in Eq. (^7|), one has 



Fi{x,Q) 



X 



2(1 - X){X, - X,)Q^ n.+Sn.=n"^^'^' 



(63) 



where in the sum all the values of constrained hy n = Hx + Uy + Uz are allowed, and ni(2) 
is defined in Eq. (^8]). By using in Eq. (^) the expressions of p'l and pi""*" obtained from 
Eqs. (^) and (|33D , J-'{n,x) is given by 



jF(n, x) 



Ipnik)- 



rvyrvz 



W[iuik)+mo] 



mo + 



Xoik'] 



kl 



u}{k) + mo 
dk 

{2Tifi'u}{kY/^ 



Xoik') 



dk 



(27r)3ecu(A;)V2 



+ 



16 



^.{k)^k^{\ - 2^ xdk') 1 



+ 



1 



ky hx 



xS') 



dk 



[uj{k) + mo] '^'^ ' (27r)3ecu(/c)i/2 ' 
It is important to note that using the completeness over rix and Uy one has 

dk 

r[n,x} = I ■ 

nx=0 n„=0 



(64) 



tpnXh) 



(27r)2 
mo + 



uj{k) + mo] 



/T[uj{k)+mo] 



XoiK. 



XoiK. 



dk. 



dk-^ 



2Tii' uj{kyl^ 
dkz 



2Txi'uo{kYl^ 



+ 



1 £' 

7==kJl - 2-) XoiK) ■ , 



+ 



ky hx 



XoiK] 



dky 



(65) 



/a'[cu(A;)+mo] '^"^ ^' 27recu(A;)V2 ' 
with Xoi^'X) = ^iK^ ^i-Y^'^'^oik'z)- Then applying similar arguments as in Appendix A, where 



mo has to be replaced by ymg + k\, one obtains, as in the pedagogical example, that only 
the negative and large values of kz contribute to the integrals in Eq. (^) in the Bjorken 
limit, and therefore ^' ^ x. Finally one can find 



- - 8a^ f dk^ , p 2 8a2p(x)(l-x) 



nx=0 ny=0 



2(27r) 



(66) 



We explicitly note that in Eq. (^) the order of lim^^^oo and of the integration over k±_ has 
been exchanged, thanks to the presence in Eq. ( |65|) of the gaussians ipoikx)ilJoiky) and to the 
power-law behaviour in k±_ of the remaining terms, for kx(y) — >■ ±cx). Equation (^) implies 
that the positive function JF(n, x) is bounded by a quantity which does not depend upon n. 
Let us come back to Eq. (^), that can be rewritten as 



Fi(x,Q) 



X 



E 

nz=M AX{0,n2—nx—ny) 



ni rii-rix 

16a2(l-x)(ni^^3^o„5o 
jF(r?, x) 



(67) 



First of all, let us demonstrate that 



lim JF(n, x) = lim JF(n, x) = 



(6J 
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Indeed in Eq. (0), there are integrals of the form 



dki ipnXki) i^o{h)gi{ki) i = x,y (69) 

where it is easy to show that 4'o{ki)gi{ki) is i) continuous, ii) bounded and iii) G £i (due to 
the gaussian form of ipo{k,j) and the power-law behaviour of gi{k.j) for ki ±00). Further- 
more also the derivative of ipo{ki)gi{ki) belongs to Ci. Therefore, we can exploit the same 
arguments presented in Appendix C and the boundedness of T{n,x) (see Eq. (|66D) to show 
that for large and Uy one has 

(J 

T{n,x) < 3,^ 3,^ (70) 

Ux Uy 

where C > does not depend upon n^. Actually, using the property of the Hermite polyno- 
mials H'^j^^{t) = 2{n + l)Hn{t) and the gaussian behaviour of ipo{ki)gi{ki), it is possible to 
show that the fall-off of JF(n, x) is faster than any power of n^Uy. 

By using the upper bound of Eq. ([70|) one can show that even for Q — 00 (i.e. 
ni(2) — >■ 00) the sum over in Eq. ( pTj) contains only a finite number of relevant terms. As 
a matter of fact, one has 

-1 ni ni—nx ni—nx—riy 

EE E Hn,x)< 

V I 2; nx=N^^i ny=0 nz=MAX{0,n2-nx-ny) 

3/2 3/2 — 

"-2; nx=Nx+l ny=0 nz=MAX{0,n2-nx-ny) ^x Uy 
ni n\—nx ^ ni oo 

E E 3/2 3/2 E -372 E -372 - ^ C^^) 

nx=Nx+l ny=0 Ux Uy nx=Nx+l f^'^ ny=0 Uy 

provided that A^^^ is large enough. By analogous arguments one can show the same property 
for the sum over Uy. Furthermore, the indices rii — rix — riy and n2 — — riy in the sum 
over can be replaced by rii and n2, respectively, by using once more the upper limit in 
Eq. (^). As a matter of fact, for large values of rii and n2 and fixed values of and Ny, 
one has 

Nx Ny "2(1) 

EE E Hn,x)< 

nx=Ony=0 riz=n2(i)-nx—ny 
Nx "2(1) ^ 

EE E 3/2 3/2 — 

nx=Ony=0 nz=n2^i)-nx-ny Ux Uy 



1 


{rii 




^2) 




1 




(ni 




^2) 




C 





( , E E ^^i^ < (72) 
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Therefore, in the Bjorken hmit, namely for large ni and n2, Eq. (|67| ) becomes 



Fi{x,Q) 



X 



V ni 



(73) 



with finite values A^^^ and Ny which do not depend upon Q. Then, adding to Eq. (|73D the 
following sum (that is vanishing for large values of A^^^ and A^j, due to Eq. ([T0|)) 



one has 



Si{x,Q) 

ni 
nz=?i2 



Fi{x,Q) 



X 



16a2(l — x){ni — 



oo oo 



E E + E E 

nx=Nx+l ny=0 nx=0 ny=Ny+l 



(74) 



X 



oo oo ni 



16a2(l -a;)(ni -na) „^on 



E E E •^(^.^) 



p(x) 



(75) 



= "2=^2 



where Eq. (^6]) has been used in the last step. 

Arguments analogous to those presented in this subsection can be applied in order to 
show that W~^^ is vanishing in the Bjorken limit also in the general case. 



5 Conclusion 

Within front-front dynamics we have studied DIS in the Bjorken limit for a relativistic 
model system composed by two spin 1/2 particles interacting via the relativistic analog of 
the harmonic oscillator potential (see Sec. 

First of all, we have explicitly calculated the DIS structure functions at ^ oo for 
our model system with the standard parton model assumptions in a Breit frame where the 
total transverse momentum P'±, as well as the transverse momentum transfer q± are zero 
(i.e., in an infinite momentum frame). We have adopted as usual a one-body electromagnetic 
current, and have shown that the structure functions are indeed given by the standard 
formulas of the parton model. It should be pointed out that a one-body choice for the 
current operator in the Breit reference frame is compatible with Poincare invariance (see 
Ref. [P). 

Then we have introduced the exact final state wave functions which properly take into 
account the FSI, and we have calculated the structure functions in the same Breit frame. In 
this frame, the components J+, J^, of the current operator have been taken in the one-body 

0, 



this 



form, while J~ has been defined through the current conservation. As shown in Ref. 
choice of the current is compatible with Poincare invariance and trivially fulfills the current 
conservation. We have shown that if one takes average values of the exact structure functions 
over small, finite intervals [xi,a;2] of the Bjorken variable x, so that {x2 — xi)/xi <^ 1, then at 
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high there exist many values of n such that Q^/ {Q^ + Sa?n) G [xi, 0:2] and in the Bjorken 
hmit such average values coincide with those given by the parton model. This averaging 
procedure corresponds to the finite resolution of the experimental measurements and allows 
to avoid some mathematical technicalities of distributions. 

It is the first time that the effects of the final state interactions in DIS have been 
exactly calculated in a relativistic model. It is worth noting that the relativistic calculation 
differs from the nonrelativistic ones considered in Ref. in several aspects. In particular, the 
Bjorken limit implies that one gets a finite contribution to the structure functions only from 
excited states with n — > 00, while the nonrelativistic approach is valid only if n <^ {mo/a)'^ 
(see Eq. (|14D). Another crucial difference between the two cases is that the nonrelativistic 
relation k^" = k'^ + qz/2 considerably differs from the corresponding relativistic expression 
that can be obtained by Eq. (34). Furthermore in the infinite momentum frame only the 
transverse components of the hadronic tensor survive in the Bjorken limit (see Sec. while 
in the nonrelativistic case the component W'^^ is the dominant one. 

One might think that the main result of this paper is a consequence of the fact that 
the states with given quantum numbers [n^, ny,nz), as well as the states with given values of 
k, form complete sets in Timt if the spin variables are dropped. Since the hadronic tensor is 
determined by a sum over all possible final states, the sum over {n^, Uy, n^) in Eq. (|^) should 
give the same result as the integration over k in Eq. (^), as a consequence of the above 
mentioned completeness. This is not the case due to the presence of the delta function in Eq. 
(|4^). Indeed, at fixed values of x and Q, the sum over {ux, Uy, Uz) is carried out only at some 
fixed values of the three indices (see Eq. (^Sf)). Then the calculation of the average values 
of the structure functions (see Eqs. ( ^Tf ) and (|59|)) involves only a small part of all values 
of 77, = 77-2. Therefore formally the completeness property cannot be used. Heuristically, one 
can say that the eigenstates of the relativistic analog of the harmonic oscillator potential 
are equivalent to the free states in the relevant part of the Hilbert space. It represents an 
interesting topic to be investigated whether the same property holds for different confining 
interactions. An analogous result on the equivalence between interacting and free eigenstates 
has been obtained for the reaction e^e~ — > hadrons in the model considered in 0. 

Our results could be considered as an argument in favor of the "common wisdom", 
according to which the FSI in the Bjorken limit is a higher twist effect (in this connection 
it could be interesting to calculate the terms of order l/Q^ in the structure functions and 
to see how they depend on the confinement radius). As noted above, our choice of the 
current is a possible choice, compatible with Poincare invariance and current conservation. 
However, these requirements do not determine the current operator uniquely and many body 
components could be present in J'^(O). Therefore one should study whether the results of 
the parton model can still be recovered in the Bjorken limit if the operator J'^(O) contains 
many-body interaction terms. This problem will be considered elsewhere within a more 
refined model, including possibly a three particle system. 
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Appendix A 

In this Appendix it will be proved that 



lim f(n, x) 



From Eqs. (|T^) and ( pU]) one has 

/(n,x) = fi{n,x) + /2(n,x) + h{n,x) + fi{n,x) + f^{n,x) 

where 

/i(n,x) = ^^(8v/^)^/V(n), 



(A.l) 

(A.2) 
(A.3) 



Un,x) = {-lT{2^f'^ <|.„(t)- 

J —CO 

. Xo(f) 2axo(x)|t|V^ 



]dt, 



(A.4) 



/3(n,x) = (-ir(2v/^)^/2/_%„(t)- 
r Xo(r) 2axoix)\t\y' 



]dt, (A.5) 



with /(ri) given by 

f{n) = / <^n{t)t'/^dt, (A.8) 



and $„(t) defined in Eq. (|T3|). 

So far no approximation is made and bi and 62 are arbitrary positive numbers. 
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In order to prove Eq. (|6lD , we have to show that for any e > there exists N such 

that 



\f{n,x) - 



^7r)i/2axo(a;) 



niQX 



3/2 



< e for all n > N 



(A.9) 



Let us first consider f^{n,x). When t = kz/a is positive and large then, as follows 
from Eq. (^21), ^ is close to 1. In turn, as follows from Eq. (^Sf), ^' — ^ is small. Therefore for 
large t we can replace ^' by ^ and then Xo(^') by XoiO^ which falls off exponentially when 
t ^ oo. Let us use the property (see Sec. 10.18 of Ref. p3| ) 



\^n{t)\<K 



(A.IO) 



where K ^ 1.086435. Therefore, for a given value of e it is possible to find 62 such that 
|/5(n, x)| < e/5 for all n. Due to Eq. (|A.10|) , the quantity 62 depends only on e, but not on 
n. 

As far as f 4^(11, x) is regarded, we note that, as follows from Eqs. (p!3D, (^HD and (|5^), 
there exist positive quantities ci{x) and C2{x) such that 



ci(x) < 



< C2(x) 



when t G [0, 62]- Now we use the fact that, uniformly in any bounded interval, 

1 r(n + l) 



$n(t) 



(2"n!)i/2 Tin/2 + 1) ^ ^ 2'^ ^ 



)] 



where = 2n + 1. Therefore, using the Stirling formula ||2^, for large n one has 



cos 



(Ar^/2^-n7r/2) + 0( 



n 



-111 



(A.ll) 



(A. 12) 



(A.13) 



Then by the Riemann-Lebesgue theorem (see, e.g., Ref. p4[) it is possible to find N\ 
such that |/4(n, a;)| < e/5 for all n > Ni. 

Let us now consider f2{n,x). When t is negative and \t\ is large, the quantity is 
small and becomes 2a\t\/mo (see Eq. (|52D ). Then, as follows from Eq. (^8|) , — a; is 
small and for t negative and \t\ large, one has 



Xo(r) 



2axQ{x)\t\''^ 



^/(^/^)l/2(^2 + ^2/^2)1/4 



niQX 



3/2 



An expansion of the difference between the above quantities as a Taylor series of 1/t^, 
shows that this difference behaves as l/|t|^/^ for large negative values of t. Therefore, using 
again Eq. ( |A.10| ), we conclude that it is possible to find hi such that |/2(n, x)| < e/5 for all 
n. The quantity hi depends only on e, but not on n. 

By analogy with the previous considerations we did for f4{n,x), it is possible to find 
^"2, such that 1/3(77., x)| < e/5 for all n > N2. 



22 



Finally, consider fi{n,x). Using Eq. (7.376) in Ref. and the Stirling formula, for 
large, odd values of n one has 



fin) = 2(-)i/^(-l)("-i)/2r(5/4)F(-(n - l)/2, 5/4, 3/2, 2) 



IT 



(A. 14) 



where F(a, b, c, z) is the hypergeometric function. In the literature (see, e.g., |2^) the asymp- 
totic behaviour of F(a, b, c, z) for \a\ — >• cx) is usually given in the case < 1. In Appendix 
B we obtain the asymptotic behaviour of F{—m, b, c, z) for c > 6 > and for z = 2. Then, 
using this asymptotic behaviour, one finds f{n) n^^'^ when n oo. 

For even values of n it is not convenient to use directly Eq. (7.376) of Ref. p5[, since 
it involves F{a, b, c, z) with b > c. However, from the relation 



i/„(t) = ^ [Hn+l{t) + 2nHn^i{t)] 



one has 



(A.15) 



fin) 



{2 



1 , ^ 



1/2 



'-Hn+i{t)^nHn^i{t) 



dt 



(A.16) 



Using again Eq. (7.376) of Ref. p5[, for odd values of the Hermite polynomial index, and 
the Stirling formula Eq. ( |A.16| ) can be cast in the following form for large, even values of n 



r(3/4)(-i)"/2 

f{n) = ^ iJ} J [(n+l)F(-n/2,3/4;3/2;2) 



[Tcn 



-r2F(-n/2 + l,3/4;3/2;2)] (A.17) 

Then by using the asymptotic expression for F{~m, b, c, 2) given in Appendix B, also 
for even values of n one finds f{n) tt^/^ when n oo. 
Therefore it is possible to find such that 

|/i(n,x)-(87r)V2^M^| <e/5 forall n > N, (A.18) 

In conclusion, if = max{Ni, N2, N^}, the condition ( |A.9| ) is satisfied and then Eq. 
(|6l| ) is proved. 

Appendix B 

In this Appendix we will investigate the asymptotic behavior of F(— m, b, c, 2) for 
m 00. In order to obtain this behavior let us use the well-known fact (see, e.g., Eqs. 
(9.111) and (8.834) of [||]) that, if c> 6 > 0, then 

F(-m, b, c, z) = I' t'-\l - t)^-'-\l - tz)^dt (B.l) 
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For z = 2, due to the factor {l — 2t)"^ only neighborhoods of t = and t = 1 contribute 
to the integral when m — >■ oo. We can replace (1 — tY~^~^ by 1 and (1 — 2t)"^ by (1 — t)^"^ 
in the first neig hborhood, and t'^'^ by 1 and (1 - 2t)'" by (_i)™t2m ^^le second one. Then 
Eq. (|B.1|) becomes 



F(-m,6,c,2) 



Tic) 



T{b)T{c-b) 



+ 



[\''-\l-t?'^dt + 
Jo 



t) 



(B.2) 



since in the first integral there is a non-zero contribution only in the neighborhood of t = 
and in the second integral only in the neighborhood of t = 1. Then, using the definition of 
the function B{x, y) (see, e.g., Eq. (8.380) of ||2^) and its expression in terms of the function 
T{z) (see, e.g., Eq. (8.384) of EB) one obtains 



r(c)r(2m + i) 
^(--'^'-'^)- r(c-5)r(6 + 2J+i) + 

r(c)r(2m + i) 

+ ^ ' r(6)r(c-6 + 2m + l) ^^-^^ 

Finally, using the Stirling formula, one finds the following asymptotic behaviour for 
m — >• oo 

F(-,n, ,,c, 2) = _ ^ ^ . (B.4) 

Appendix C 

In this Appendix it will be proved that W'^'^ is vanishing in the Bjorken limit. 
Since from Eqs. (ISO) and (|5B|) one has = C,'^Q{^ ~ x)/x and from the delta 



function on the masses (see, e.g., Eq. (|55D ) one obtains Q = 2ay2n^x/{l — x), then apart 
from constant factors the relevant integral in the calculation of W^^ is 



e (t2 + m2/a2)i/4 

/oo rb 
<l>n{t)^{t)dt = Vn lim / <i>n{t)^{t)dt (C.l) 



with ^(t) = ^i^/^') Xo{i')/{t^ + m2/a2)i/4. Let us investigate the limit 

lim / $„(t)^(t)dt (C.2) 

n— >oo J— c 
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We will show that this limit is uniform with respect to the extrema. First of all, let us observe 
that the function \Ef(t) is i) continuos, ii) bounded and iii) G £i. Indeed, it falls exponentially 
for t ^ oo and as for t — > — oo (as it can be argued from the discussion in Appendix 

A). Furthermore also the derivative of \E'(t) G Ci, due to the behavior of \E'(t) for t —>■ ±00. 
Since we have a bounded interval in ( |C.2|) , we can use the asymptotic expression ( |A.13| )) for 

= i—y/'[cosiN'/H - ^) + 0(n-V2)] 
nn I 



with = 2n + 1. The second term in Eq. (|C.3D gives an integral in Eq. (|C.2|) that uniformly 



vanishes as with respect to the integration extrema, due to the property iii). Also the 

first term in Eq. ( |C.3| ) produces an integral in Eq. ( |C.2| ) that uniformly vanishes as l/n^/"^ 
with respect to the integration extrema. This can be shown with an integration by parts 
and exploiting the property ii) and the integrability of \d'^(t)l dt\. 

Therefore the limit in Eq. ( |C.2| ) is zero and it is uniform with respect to the integration 



extrema. Then one has 

lim lim / <l>„(t)^(t)rft = 

c— >oo 

lim lim ^ t $„(t)^(t)rft (C.4) 

n^oo h—*oo J 
c^oo 

and finally 

lim 4 = (C.5) 



n— >oo 
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